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This article presents the results of a study of a class of error-correcting codes 
called partial-unit-memory convolutional codes, or PUM codes for short. This class 
of codes, though not entirely new, has until now remained relatively unexplored. 
This article shows that it is possible to use the well-developed theory of block codes 
to construct a large family of promising PUM codes. Indeed, at the end of the 
article the performances of several specific PUM codes are compared with that of 
the Voyager standard (2, 1,6) convolutional code. It was found that these codes can 
outperform the Voyager code with little or no increase in decoder complexity. This 
suggests that there may very well be PUM codes that can be used for deep-space 
telemetry that offer both increased performance and decreased implementational 
complexity over current coding systems. 


I. Introduction 

This article gives a general construction for, and sev- 
eral interesting examples of, partial-unit-memory (PUM) 
convolutional codes. First, some definitions and notation 
are established. 

A convolutional code C of length n and dimension k 
over a field F is defined by an encoder 

G{D) — Go + G\D + ■ * • d* G\fD (1) 

3 Consultant to the Communications Systems Research Section from 
the California Insitute of Technology. 

2 Independent consultant to the Communications Systems Research 
Section. 


where Go, Gi , . . . , G M are k x n matrices with entries 
from F 3 The ratio R = kjn is called the rate of the 
code, and M is the code’s memory. If u(D) = «o + «i D 
+ u 2 D 2 + • • • is the input to the encoder (where the uP s 
are elements of F> and D is an indeterminate), then 
X (D) - u(D)G(D) is the output, which is also called a 
codeword. The encoder is said to be noncat astrophic if 
no infinite-weight input produces a finite- weight output. 
The free distance, c?f ree , o{ a convolutional code C is de- 
fined to be the minimum weight of any nonzero codeword 
x(D) E C. If the encoder G{D) is noncatastrophic, then 
df r ee is the minimum weight of all codewords u(D)G(D) 
generated by inputs u(D) of finite weight. All other things 
being equal, it is generally desirable to have the quantity 

3 In tliis article, it is always assumed that F = GF{ 2), but most or 
all of the results generalize easily to other finite fields. 
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Q = Rdf ree as large as possible, since Q is the asymptotic 
coding gain of the code, which is a good measure of the 
communications improvement afforded by the use of the 
code. 

A convolutional code C has state complexity m if the 
sum of the maximum degrees of the rows of G(D) is m, 
This terminology reflects the fact that a physical encoder 
for C based on G(D ) has 2 m states. It is desirable to have 
m as small as possible, since the computational complexity 
of the Viterbi decoding algorithm for C is proportional to 
2 m . 

The notation “[n, k, m ) d\ code” is introduced to de- 
scribe a convolutional code of length n, dimension 
state complexity m, and free distance d. The notation 
“(n,&,m)” code is sometimes used to describe the same 
code without explicitly referring to its free distance. For 
example, in this notation, an [n, k , d] block code, which can 
be viewed as a convolutional code with M = 0, is both an 
[n,fc,0,d] and an (n,fc,0) convolutional code. For a given 
n, k } and m, a code for which c/f ree is as large as possible 
is said to be an optimal (more properly, distance optimal) 
code. 

A convolutional code with M = 0 is just a block 
code. A convolutional code with M = 1 is called a unit- 
memory convolutional code. Unit-memory codes seem to 
form a class that lies halfway between block and con- 
volutional codes. They were first studied seriously by 
Lee [5], who found a number of interesting examples of 
unit-memory convolutional codes. Thommesen and Juste- 
sen [10] have obtained bounds on the performance of unit- 
memory codes, and Justesen, Paaske, and Ballan [3] have 
constructed a class of unit-memory codes which they call 
quasi-cyclic codes. This article studies another subclass of 
unit-memory codes called partial-unit-memory codes. 

For a unit-memory convolutional code, the state com- 
plexity is just the number of nonzero rows of G\ . If some 
of the rows of G i are zero, i.e., if m < £, then it is said 
that the code is a partial-unit-memory (PUM) convolu- 
tional code. Partial-unit-memory codes were introduced 
by Lauer [4], who constructed several optimal PUM codes. 
Some general constructions and further examples of PUM 
codes (under the name finite-state codes) were given in [7] 
and [8]. This article should be viewed as a continuation 
of these earlier studies, in which, among other things, it 
is shown that many of these earlier results follow from the 
authors’ methods. (For example, Lauer’s equidistant PUM 
codes appear in Example 2 in Section III of this article, 
and the general construction of Theorem 5 in [7] appears 
as Corollary 4 in Section II.) 


Here is a summary of this article. Section II gives a 
general construction for PUM codes based on the exis- 
tence of certain block codes. Roughly speaking, the main 
result is that if there exist two distinct [n,k,d 0 ] block 
codes with a common [n,fc*,d*] subcode, then there ex- 
ists a noncatastrophic [n,Jfc,ifc — Jb* , d] PUM code with 
d > min(d*,2d 0 )« (As a point of comparison, Theo- 
rem 5 in [7] shows that if there exists a single [n,jb,d 0 ] 
block code with an [n,^*,cT] subcode, then there exists 
a noncatastrophic [n, k - 1, k - jfc* ~ 1, d\ PUM code with 
d > min(<i*, 2do).) Since there is a huge existing cata- 
log of block codes, what this means is that it is possi- 
ble to construct a very large number of interesting PUM 
codes. This is illustrated in Section III with several ex- 
amples called Hamming, Reed-Muller, and Golay PUM 
codes. While these examples are possibly interesting and 
potentially important, the authors believe that they have 
only scratched the surface, and hope that future authors, 
using these techniques, or ones of their own devising, will 
unearth many more examples. 


II. Main Results 

Theorem 1 . Suppose that C 0 is an [n y k t d 0 ] lin- 
ear block code, C\ is an [n } k,di] linear block code, and 
Co 7 ^ Ci . Suppose further that Cq and C\ contain a common 
subcode C*, which is an code. Then there ex- 

ists a noncatastrophic [n, k, m, d\ PUM convolutional code, 
with m = k - k* and d > min(d*, do -f c/i). 

Proof; Begin by choosing k x n generator matrices G 0 
and G\ for Cq and C\ of the form 



where K* is a k* x n generator matrix for C* . Note that 
both K 0 and K x are m x n matrices; for future reference, 
let Cq and C\ be the corresponding codes, i.e., 


o 

II 

*© 

(3) 

£ 

II 

* _ 

(4) 

Next, define the matrix Gq as 


<*-[£] 

(5) 
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where O is a k* x n matrix of 0’s. Then define a k x n 
polynomial matrix G(D) as follows: 


which means that x(D ) is a nonzero vector in the rows- 
pace of K*> i.e., a nonzero codeword in C*, and so weight 
(x(D)) > d * in this case. 


G(D) = Go + G\D (6) 


Plainly, G(D ) is the generator matrix for an (n, k , m) PUM 
code. The proof will be complete when the following 
two things are shown: (1) d( ree > min(d*,do + an< i 
(2) G(D ) is noncatastrophic. Begin with the assertion 
about dfree- 

Assume that u(D) = uo + u x D + u^D 2 + * • * is a finite 
nonzero input sequence, where uo / 0 a nd each u t - is a 
fc-dimensional row vector. Then the corresponding output 
sequence is x(D) = u(D)G(D) — xo + x x D + x^D 2 + * * ’» 
where 


Xo = UqGo 

Xi = UiGo + for i > 1 


(?) 


If ti = (^1,^2, • • • I L*k) is a ^-dimensional vector, u L (the 
left part of u) and u R (the right part of u) are defined as 
follows: 


If, on the other hand, / 0 for some i > 1, 

let M denote the largest such index. Then, u = 0, and 
Eq. (9) implies 

x(D) = uqGo + * ■ * + u m- \\G\D M (12) 

But u 0 Go is a nonzero word from Co, and so has weight 
> d 0} and is a nonzero word from C Xi and 

so has weight > d x \ thus, weight (x(D)) > do 4* d x in this 
case, which proves the assertion about d( ree . 

It remains to be shown that G(D ) can be chosen non- 
catastrophically. Lemma 1, which follows, tells, in princi- 
ple, whether a given G(D) is catastrophic or not. Lemma 2 
then tells that it is always possible to choose the matrices 
Go and Gi so that G(D) is noncatastrophic. 

Lemma 1. Let the linear transformation T : Co — * ► C x 
be defined by itGo — ► uGj. Then G(G) is noncatastrophic 
if and only if every subspace of Co fixed by T is a subspace 
of C*. 


u L - (/ii, . • ■ 

u R = . ,/ifc) 


Proof: Denote the rows of K * by (x x , Xj, . . . , x*.), the 
rows of K o by (y x , yi, • • • , 3/m)> and the rows of K x by 
( 8 ) (z 1 ,z2,...,Xm)- Then T is completely characterized by 

the k values 


Then [see Eq. (2)J f for any vector u, one has 


uG 0 = u l K* + u r K 0 , uG? = u r K i (9) 

After combining Eq. (8) with Eq. (6), one has 

Xo = Uo K* + u R Ko "1 

( 10 ) 

x i = u?K 0 + [ut,u?_ l }G u fori>l J 

Now, either = 0 for all i > 1, or not. It will be 

shown that weight (x(D)) > d * in the first case, and that 
weight (x(D)) > d 0 + d x in the second case. 


Tx i = x i, 

for i — 1,2,.. 


(13) 

Tyi = Zi, 

for i — 1,2,.. 


(14) 


Note that Eq. (13) says that T not only fixes C\ it fixes 
C* pointwise. 

It is first assumed that every T-fixed subspace of Co 
is a subspace of C*, and then shown that G(D ) is non- 
catastrophic. Let u(D) be a nonzero input such that the 
corresponding output x(D) is finite, i.e., x* = 0 for i > 

If one defines 


= ec* 


(15) 


If = 0 for all i > 1, then by Eq. (6) and 

Eq. (9), one has 


bi - u r Kq E Cq (16) 


x(D) — uqGq = Uq K* 


(11) one has, by Eq. (9), 
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Xi = di -f 6 * + TVi, for i > 1 (17) 

Thus, since it is assumed that = 0 for i > to, it follows 
that (recall that the codes are binary) 

Tbi-i = bi + cLi y for i > * 0 (18) 

so that (6i 0 ,6i 0 +i, . . . ) + C* is a T-fixed subspace of Cq. 
However, it is assumed that all T-fixed subspaces of Co are 
subspaces of C* , and so bi E C * for all i > i 0 . But since 
the rows of Kq are linearly independent of the rows of K * , 
this means that &,* = 0, and so uf- — 0, for all i > t 0 . But 
then, by Eq. (17), a,- = 0, and so uf = 0, for all i > i 0 . 
Thus, the input u(D) is necessarily finite, and so G(D) is 
noncatastrophic. 

Conversely, suppose that B is a nonzero T-fixed sub- 
space of Co that properly contains C*. Choose a q E C* 
and b 0 E B — Cq arbitrarily. Now, since B is T-fixed, Tb 0 
is also in £?, and so it can be decomposed uniquely into 
the sum of an element of C* t which is called ai, and an 
element of Cq, which is called b\. Note that since b 0 ^ 0, 
then bi ^ 0 also, for otherwise T would map the k* + 1 
dimensional space C* + ( 6 0 ) into the k * -dimensional space 
C*. This process is continued inductively by constructing 
an infinite sequence of pairs (a t ‘,6j), with a, E C* , bi E Cq, 
&i 7 ^ 0 , such that 


Tbi = a,- + i + 6 »+i, for i > 0 (19) 

Now for each i > 1, define the vector u, as follows: 


ufK* = a, 

(20) 

uf K 0 = ^ 

(21) 

Since by / 0, then by Eq. (20), ^ 0, and so the sequence 

(uy) is an infinite sequence of nonzero elements. It will now 
be shown that, if (ti,) is the input, then the corresponding 
output is finite. One has 

j? 

n 

9 

-S3 

g 


= r([o,«f]G 0 ) 


= [ 0 ,uf]Gi, 

by definition of T 


= of Ki 


and so by Eq. (9), for i > 1, 

Xi = ufK* + tif + of.,*!, by Eq. (9) 

= Qi + + Tbi^i, by Eqs. (19) and (20) 

= by Eq. (18) 

Thus, the infinite input sequence (u;) produces a finite 
output sequence (x,-), and so G(D) is catastrophic, as was 
asserted. □ 

Corollary 1 , If Co H Ci = C% then any generator 
matrix of the form Eq. ( 6 ) is noncatastrophic. 

Proof: If B is a T-fixed subspace of T, then since 
T : Co — ► Ci and B = T(H), B is a subspace of both Co 
and C\. Thus, B C Co D C\ ~ C*, and so by Lemma 1 , 
G(D ) is noncatastrophic. O 

Lemma 1 allows one to tell, in principle, whether or 
not a given G(D) is catastrophic. Corollary 1 assures that 
if Co fl Ci = C* y then nothing can go wrong. However, if 
Co H Ci DC*, more work is necessary to find a noncatas- 
trophic generator matrix. Lemma 2, which follows, gives 
an explicit construction for a noncatastrophic G{D) in the 
general case. 

Lemma 2. Suppose that Co and C\ are subspaces 
of V n (F) } the n-dimensional vector space over F } with 
Co ^ Ci but dim(Co) = dim(Ci) = k t and that C* is a 
subspace of both C 0 and C i} with dim(C*) = k*. Then, if 
(«i , , tLk * ) is a basis for C* , there exist bases for Co 

and Ci of the form 


(Cq) (tii , . . . f ujfc* , 

• > a m) 1 

( 22 ) 

(Ci) = («!,...,«*•, 

\pm) ) 


such that the linear transformation T 

' Co — ► Ci 

defined 


by 

Tui = ti,*, for * = 1 ,... } k* } 

> (23) 

Toti - Piy for ? = l,...,m J 

fixes no subspace of Co larger than C*. 

Proof: Begin by constructing two descending se- 

quences of subspaces (.4*) and (Bi): 
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Co = Aq D A\ D-O -4jv+i — £ 
Ci = B 0 DBiD--OBjv+i = C* 


such that 


the vectors (tii , . . . , ujfc* , ai , . . . , cx m ), and let the rows of 
G\ be the vectors (ui, . . . , u** , /3i ( --*iAn) i n Eq. (21). 
Then, the mapping T : Co — > ► C\ defined by uGq — * uG\ 
is the same as the mapping described in Lemma 2, and so 
the resulting G(D) is noncatastrophic. This completes the 
proof of Theorem 1. ^ 


dim(v4i) = dim(yBj), i = 0,1,. . . , AT + 1 

Ai + 1 = Ai n Bi } i = 0> I* • • • * N 

Figure 1 illustrates the construction of Lemma 2. This 
construction can be done inductively as follows. Assume 
that Aq 9 Ax, ...,Ai and BQ,Bu...,Bi have already been 
constructed. (For i = 0, this simply requires setting Ao = 
Co and Bq = C\.) Let Ai+i = A{C\B{. If Ai+\ = C , define 
Bt+i = C*, N = i, and stop. Otherwise, one has C* C 
w4*+i C Bi, and so by Lemma A2 in Appendix A, there 
exists a subspace Bi + 1 ^ ^4i+i such that C C Bi+ \ C $*, 
with dim5i+i = dimA'+i- 

Now define integers fro? • * ■ j by ki = dim(*4i) — fc , 
so that 

< fc/v-i < • ■ • < — m 

Next, using Lemma Al in Appendix A, choose bases 

for C 0 and (ui, . .. , iff , 

/? m ) for Ci such that 


{« 1,. • 

«1,-' 

-,OCk,) = Ai 


■,Uf, Pi,-- 

II 


for i — l,2,...,iV\ Now define the transformation T as 
in Eq. (22). Plainly, T fixes C* pointwise, and also, from 
Eq. (23), 

T:Ai-+B iy for 1 = 0,1,..., AT (25) 

If now £> is a subspace of Co fixed by T, then V C Co = 
>to) and = T('P) C T(Co) = Cx = Bq, so that P C *4o D 
Also, P = T(P) C T(Ai) = B u using Eq. (24), 
so that P C Ai O Bi = A 2 . Continuing inductively, one 
finds that in fact P C ^ 43 , . . . , P C 1 = C*. Thus, 
a linear transformation T is constructed such that any T- 
fixed subspace of Co is a subspace of C*. O 

Lemma 2 tells how to construct noncatastrophic gen- 
erator matrices Go and Gx: Just let the rows of Go be 


Corollary 2. Suppose that Co is an [n,fc,do] linear 
block code, and that C* is an [n,Jr*,<P] code, which is a 
subcode of C 0 . If the automorphism group of C* contains 
a permutation that does not fix Co, then there exists an 
[n, Jb, m,cf] PUM convolutional code, with m — k — k* and 
d > min(d*,2do)- 

Proof: Let tt be an automorphism of C* that does not 
fix C 0 , and let Ci = CJ. Then, C x is an [n,Mo] code not 
equal to Co- Now apply Theorem 1. O 

Corollary 3. If C 0 is an [n, Jk, d 0 ] linear block code that 
contains the all-ones vector, and if fc ^ l,n — l,n, then 
there exists an (n,fc,fc — 1) PUM code with d( ree > 2do. 

Proof: Here, Corollary 2 is applied, with C* being 
the [n , 1 , n] code consisting of the two vectors [00 • • •] 
and [11*1]. Clearly C* is fixed by all permutations of 
{1,2,..., n) - Furthermore, the only binary linear codes 
that are fixed by all permutations of {1,2, . . . , n) have di- 
mensions 0, 1, n - 1, or n, so there must be an automor- 
phism of C* that doesn’t fix.Co- Thus, by Corollary 2, there 
exists an (n,fc,Ar— 1) PUM code with df re e ^ min(2do, n). 
However, since k > 2, the mimimum distance do of Co must 
be < n; and since Co contains the all-ones vector, there 
must be a word of weight n — do. Hence, n — d 0 < do, and 
so do < n/ 2. Hence, min(2do» n ) = 2do, so that in fact 

> 2<f 0 - a 

Corollary 4. (Same as Theorem 5 in [7]). Suppose 
that Co is an [n,Jb,d 0 ] linear block code, C* is an [n,fc*,cT] 
code that is a subcode of Co, and k — k* 2. Then, for 
every integer i in the range 1 < t < k — k* — 1, there is 
a noncatastrophic [n, k — t, k — i — k* , d\ PUM code with 
6 > min(d* , 2 d 0 ). 

Proof: Let Cq be any (n, k — i) subcode of Co that con- 
tains C*. (The conditions on i guarantee that dimC*z < 
dim Cq < dim C 0 , so this is possible.) By Lemma A2, there 
exists a subcode C[ not equal to Cq but having the same 
dimension, and also lying between C 0 and C*. Thus, Cq 
and C[ are both [n.k - i y d'] block codes, with d' > d 0 . 
By applying Theorem l’to the codes Cq, CJ, and C , one 
obtains a noncatastrophic [n, k - z, k - i - k * , d] PUM code 
with <5 > min(d",2d') > min(</*,2rf 0 ). O 


61 


III. Examples 

In this section, four examples of PUM codes are pre- 
sented, that were constructed with ’the help of the results 
in Section II. Example 1 describes a Hamming [8, 4, 3, 8] 
PUM code, which was originally discovered by Lauer [4]. 
Example 2 gives a generalization of Example 1 to a class 
of Reed-Muller [v2 + 1,/i, j/2* 1 ] PUM codes, one code 
for each pair of positive integers (/i,^) except (1,1) and 

(2. 1) . (The code of Example 1 corresponds to the pair 

(3. 1) .) The codes in Example 2 were also found, using 
different methods, by Lauer. Finally, in Examples 3 and 4, 
two new PUM Golay codes, with parameters [24, 12, 7, 12] 
and [24, 12, 10, 16], are presented. 

Example 1 (a Hamming PUM Code). Let Cq be 
the [7,4,3] binary cyclic code with generator polynomial 
go(x) = 1 -f x + x 3 , and let C[ be the [7, 4, 3] binary cyclic 
code with generator polynomial g\{x) = l + x 2 + x 3 , Take 
as a generator matrix for Cq the 4x7 binary matrix Gq, 
whose rows are g l (x)g 0 (x), go(x), xg 0 (x) } and x 2 g 0 (x) t 
and for C[ the 4x7 binary matrix G\ , whose rows are 
ffo(*)l7i(ar), 0i(*), xg x (x), and x 2 gi(x ), i.e., 

( 1 1 1 1 1 1 1\ 

110 10 0 0 
0 110 10 0 
0 0 110 10/ 


( 1 1 1 1 1 1 1\ 

10 110 0 0 
0 10 110 0 
0 0 10 110/ 

Now, if each code is extended to length 8 by appending 
an overall parity-check, one obtains codes Cq and C\ } 
both of which are binary [8,4,4] codes with generator 
matrices 

( 1 1 1 1 1 1 1 1\ 

1110 10 0 0 
10 110 10 0 
10011010/ 


/I 1 1 1 1 1 1 1\ 

r 1 1 0 1 1 0 0 0 
1 110101100 
\1 0 0 1 0 1 1 0 / 

Since C' 0 D C[ is the binary [7,1,7] repetition code, it 
follows that CoflCi is the binary [8,1,8] repetition code. 
Thus, in Theorem 1, C * can be taken to be the [8,1,8] 
repetition code with the 1x8 generator matrix 

IC = ( 1 1 1 1 1 1 1 1 ) 

It follows from Corollary 1 that the matrix 


G(D) 


1 1 1 1 1 1 i 1\ 

l + D 1 + D 1 D 1 + D 0 0 0 

1 + D 0 1 + D 1 D l + D 0 0 

\1 + D 0 0 l + D 1 D l + D 0/ 


generates a noncatastrophic [8, 4, 3, 8] PUM code. Further- 
more, from [4] (Formula (3) with L = 0) or [7] (Corollary 1 
to Theorem 1, with L = 1), any (8, 4, 3) convolutional code 
must have cff ree < 8, so this code is optimal. 4 □ 

Example 2 (Some Reed-Muller PUM Codes). 
Let /i and v > 1 be positive integers. Let A ^ be the 
[2^,/i-h 1,2^ _1 ] first-order Reed-Muller code, and let B p 


4 Tliis code first appeared in the literature in [4], Table 1. It is ap- 

parently used by the Soviets in their Regatta space communication 
system. 


be the [2^, 1,2 M ] zeroth-order Reed-Muller code (a repeti- 
tion code), which is a subcode of A ^. 5 Now let Co(/i,^) 
be the [v 2^,/i-f l,t/2 M_1 ] code obtained by repeating A ^ 
v times, and let U*(/i,z^) be the [t^, 1, v 2^] code obtained 
by repeating B M v times. Then, according to Corollary 3, 
unless = (1,1) or (2,1), there exists a noncatas- 

trophic [u2 ^ , /i -f l,/i,i/2^] PUM code. These codes are 
all optimal by the above-cited bounds in [4] or [7]. (This 
family of codes was originally constructed by Lauer [4], us- 
ing a different approach. He called them equidistant PUM 


5 Mac Williams and Sloane [6], Chapter 13, is a good reference for 
Reed-Muller codes. 
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codes. A similar family of [2",/x,A* - 1,2*] codes was con- 
structed in [7], Example 4.) ^ 

Example 3 (the [24,12,7, 12] Golay PUM Code). 

It is well known that there exists a [24, 12,8] binary linear 
code; viz., the famous Golay code. It turns out that there 
are two isomorphic copies of the Golay code that contain 
a common [24, 5, 12] subcode, so that by Theorem 1, there 
exists a noncatastrophic [24, 12, 7, 12] PXJM convolutional 
code. In this example, the construction of this code is 
detailed. 

Define, for A,B,C,D,E,F <E GF{ 8), the following two 
functions: 

y ) = Tr(Axy) + Tr((B + Cy)* 6 ) (26) 

9 d,e,f{ x , V ) = Tr((-Oy + + Tr((Fy)x 6 ) (27) 

where Tr(x) = x+x 2 + x 4 is the trace mapping from GF{ 8) 
to GF(2). Then, if /? is a fixed nonzero element of trace 
0 in GF( 8), the following set of 2 12 length-24 vectors is a 
[24, 12,8] Golay code, which is called Aq: 

Ao = \fA,B,c{ x ,P) + Col fA,B,c{x,P 2 ) + e l I 

x fA,B,c{ x >P*) + f 2] reGF (g) ( 28 ) 

In Eq. (28), the parameters A, B, and C assume all values 
in GF(8), and the parameters to, Gi and C2 assume all 
values in GF(2). The proof that A 0 is indeed a [24, 12, 8] 
code appears in Appendix B as Lemma B5. 

Similarly, the following set of 2 12 length-24 vectors is 
another [24, 12,8] Golay code, which is called B 0 : 

B 0 = [gD,E,F( x >P) + 6 o\9d,e,f( x ,0 2 ) 

+ ?>i\gD t E,F( x i p 4 ) + (29) 

In Eq. (29), the parameters D, E, and F assume all values 
in G F(8), and the parameters So, Si, and S 2 assume all 
values in GF(2). The proof that B 0 is indeed a [24, 12, 8] 
code appears in Appendix B as Lemma B6. 

In Appendix B (Lemma BIO), it is shown that Ao Cl Bo 
is a [24, 9, 8] code consisting of the following set of vectors: 


A, = [fA,o,c( x ,P) + eolW(*,/? 2 ) + f il 

X /^,0,C(Z.^ 4 ) + ^ 2 ] x £<7F(8) ( 30 ) 

The code Ai, in turn, contains a [24,5, 12] subcode con- 
sisting of the following set of vectors: 

A3 = [/o,o,c(z./?) + fo|/o,o,c(z,/? 2 ) + fi | 

X /o, 0 ,c(z,/? 4 ) + Co + Cl] ieGF (g) (31) 

(see Lemma B8 in Appendix B). Finally, A3 contains a 
[24, 2, 16] subcode A 4 : 

A 4 = [fo|ci|co + Cl] r gGf’(8) (32) 

(see Lemma B9). It follows from Theorem 1 that there 
exists both a [24, 12,7, 12] code and a [24, 12, 10, 16] code, 
and by the bounds in [4] and [7], they are optimal. 6 To 
actually construct noncatastrophic generator matrices for 
these two codes, however, more work is necessary. Here 
are the needed intermediate subspaces (see Fig. 2): 

Bi = [go,E,F( x ,0) + 6 o \go,EA x >0 2 ) + M 

x go,E,F{ x ,P 4 ) + 5 2 ] reGF( g) (33) 


A2 = [/o,0,c(z,/?) + £o|/o,0,c(z,/? 2 ) + fll 

x /o,o,c(z,£ 4 ) + c 2 ] x€CF(8) (34) 


Bi = [go,efi{ x ,0) + So\go,e,o( x ,P 2 ) + <5l| 
x yo,e,o(z,/? 4 ) + ^2] x€ OF(S) 

In Eq. (35), e assumes only the two values 0 and 1. 


6 In [8], [24,5,12] and [24,2,16] subcodes of a [24,12,8] Golay code 
were found, which led, via Corollary 4, to the construction of both 
[24,11,6,12] and [24,11,9,16] PUM codes. 
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In order to show that the subspaces in Fig. 2 behave as 
depicted, the following must be proved: 

AqC\ Bo = A\ (36) 


A l nB 1 = A 2 (37) 

a 2 n B 2 = A% (38) 

These relationships are proved in Appendix B in Lem- 
mas BIO, Bll, and BI2. 

It thus follows that for the [24,12,7,12] code, a non- 
catastrophic choice for Go and G\ is as follows: 


’ /ooi 


" /ooi 

hop 


foop 

/oO/P 


loop 2 

Oil 


Oil 

101 


101 

100 

/ioo 

Gi = 

#010 

100 

o 

o 


9opo 

//POO 


<7o/Po 

/oio 


9 ioo 

fopo 


9 poo 

— t 
o 

c* 


-£/poo- 


Here, fA,B,c denotes the length-24 vector obtained by tak- 
ing e 0 = ci = c 2 = 0 in Eq. (28), and gA,B,c denotes the 
length-24 vector obtained by taking <5 0 = = <52 = 0 in 

Eq. (29). The first five rows of G 0 and G\ are identical, 
and they generate the [24,5, 12] subcode referred to above. 
In binary, these two matrices are as follows: 


'01110100 

00111010 

01001110 ' 

00111010 

10011100 

10100110 

10011100 

01001110 

11010010 

00000000 

11111111 

11111111 

11111111 

00000000 

11111111 

11111111 

00000000 

00000000 

00101110 

01011100 

01110010 

01011100 

10111000 

11100100 

10111000 

01110010 

11001010 

11101000 

11101000 

11101000 

01110100 

01110100 

01110100 

.00111010 

00111010 

00111010 . 


'01110100 

00111010 

01001110 ' 

00111010 

10011100 

10100110 

10011100 

01001110 

11010010 

00000000 

11111111 

11111111 

11111111 

00000000 

11111111 

10010110 

10010110 

10010110 

11111111 

00000000 

00000000 

00101110 

00101110 

00101110 

01011100 

01011100 

01011100 

00101110 

01011100 

01110010 

01011100 

10111000 

11100100 

.10111000 

01110010 

11001010 . 


□ 

Example 4 (the [24,12,10,16] Golay PUM 
Code). The [24,5,12] binary linear code of Example 3 
contains a [24, 2, 16] subcode, so that by Theorem 1, there 
exists a noncatastrophic [24, 12, 10, 16] PUM convolutional 
code. In this example, the construction of this code is de- 
tailed. The subspaces A 0 , B 0} A x , B u A 2 , and A 4 defined 
in Example 3 are used. Additionally, subspaces B' 2) A 3) 
and B 3 are defined as follows: 


B*2 — + <M0o,£.o(z j /? 2 ) + £>i\9q,e,o{x,P a ) 

(39) 


A 3 — fcoklk2]reGF(S) (40) 


B's - [<7o,«,o (*,/?) T S Q \g 0je 0 (x f /3 2 ) 4- 6\ too.e.oO*, P 4 ) 

+ ^0 + «i] reG/ . (8) (41) 


In Eq. (41), e assumes only the two values 0 and 1. 

The proof that the subspaces behave as depicted in 
Fig. 3, i.e., that A 2 fl B ' 2 = A 3 and A ' 3 n B f 3 = A 4) is 
given in Appendix B, Lemmas B13 and B14. 

Now one can see that a noncatastrophic choice for Go 
and G\ for this code is as follows: 
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Oil 

101 

100 

/ 001 
foop 
foop 3 

f 100 
fp 00 
fp 3 00 
/ 010 
fopo 
fop 3 o 


Oil 

101 

g oio 

100 

9opo 

9op 3 o 

f 001 
foop 
foop 3 

f 100 
fp 00 
fp 3 oo 




In binary, these matrices are 

'00000000 11111111 11111111' 
11111111 00000000 11111111 
11111111 00000000 00000000 
01110100 00111010 01001110 
00111010 10011100 10100110 
10011100 01001110 11010010 
Go ~ 00101110 01011100 01110010 
01011100 10111000 11100100 
10111000 01110010 11001010 
11101000 11101000 11101000 
01110100 01110100 01110100 
.00111010 00111010 00111010 . 


'00000000 11111111 11111111 " 

11111111 00000000 11111111 
10010110 10010110 10010110 
11111111 00000000 00000000 
00101110 00101110 00101110 
01011100 01011100 01011100 
Gl = 01110100 00111010 01001110 
00111010 10011100 10100110 
10011100 01001110 11010010 
00101110 01011100 01110010 
01011100 10111000 11100100 
.10111000 01110010 11001010. 

□ 

The codes in Examples 1, 3, and 4 are quite interest- 
ing as combinatorial objects, but they have potential for 
applications. To illustrate, Fig. 4 shows a plot of the per- 
formance of these three codes and the NASA standard 
[2, 1,6, 10] (non-PUM) code on an additive white Gaussian 
channel. Figure 4 shows that the low-complexity Ham- 
ming [8, 4, 3, 8] code is only a bit weaker than the NASA 
code, while the two Golay codes are both a bit stronger. 
Since the state complexity of the [24, 12, 7, 12] Golay code 
is only 1 greater than that of the NASA code, it may be 
that there is a relatively low-complexity decoding algo- 
rithm for this code, whose performance will significantly 
exceed that of the NASA code. In any case, these perfor- 
mance curves certainly justify a serious study of efficient 
decoding algorithms for these and other PUM codes. 
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1 = C* 


Fig, 1. The construction of Lemma 2. 


Fig. 3. The subspaces needed for the construction 
of a noncatastrophlc encoder for the [24, 12, 10, 16] 
code. 



Fig. 2. The subspaces needed for the 
construction of a noncatastrophlc en- 
coder for the [24, 12, 7, 12] PUM code. 


□ [2, 1,6, 10] (NON-PUM) 
• [8, 4,3,8] 

■ [24, 12, 7, 12] 

O [24, 12, 10, 16] 


E b !N 0 , dB 

Fig. 4. Performance curves for three PUM codes, compared with 
the NASA standard [2,1,6,10] code, on an additive white 
Gaussian channel. 




Appendix A 

Two Results From Linear Algebra 


In this Appendix, two simple results from linear algebra 
are provided that are needed in the proof of Lemma 2, 

Lemma Al. If (0) = Ao C A\ C C A m — V 
is an ascending chain of subspaces of an n-dimensional 
vector space V, with dim *4** = ft*, then there exists a basis 
(a!,..., a n ) for V such that 

(ai for i= 1, — (A-l) 

Proof: One proceeds recursively, as follows. Choose a 
basis , . . . , a^) for A\* If m = 1, one is done. Oth- 
erwise, by using a standard result in linear algebra [ 2 , 
Lemma 4.2.5], the basis {<* 1 , . . . , OfjbJ for A\ can be ex- 
tended to a basis (<* 1 , . . . , ot , . . . , »Jb a ) of A 2 , etc. d 

Lemma A2. Suppose that V is an n-dimensional vec- 
tor space over F, and S and T are subspaces of V with 
S C T C V. Then there exists a subspace 7' ^ T such 
that dimT' = dimT and S C T* C V. 

Proof: Suppose that dim S = k and dimT = k + j , 
where j > 0. By Lemma Al, it is possible to find a basis 
for V of the form 


where k -f j + ft = n, and 

(<*!,... ,C*jfe) = S 
(ai, . . . .orjbj/Ji, . . - ^ 

If ft > j, define T' as follows: 

T' = 

If, on the other hand, ft < j, define T f as follows: 

T' = <cvi , . . . , a* , 71 , . . . , 7h. £1 ? • • • > fy - * ) 

□ 
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Appendix B 

Proofs Needed in Examples 3 and 4 


In this Appendix, the assertions made in Section III 
about the subspaces A 0 , A u A 2 , A 3 , A 3 , A 4 , B 0 , B u B 2} 
B f 2) and B 3 of V-m( 2) are proved. 

Lemma Bl. Let A and B be elements of GF{ 8 ) such 
that Tr(Ax + Bx 6 ) = 0 for all x £ GF( 8 ). Then A — B = 
0 . 

Proof: Since Tr (y) - y -f y 2 -f y 4 and y 8 = y for all 
y £ GF(8) y it follows that 

Tr(Ax -f Bx 6 ) = Ax -f A 2 x 2 + A 4 x 4 + Bx 6 -f B 2 x 5 H- B 4 x 3 

(B-l) 

for all x £ GF(8). Thus, the equation Tr(Ax 4- Bx 6 ) = 0 
is a polynomial equation of sixth degree with 8 roots in 
GF(8) } and hence, the coefficients of the polynomial must 
be zero, i.e., A — B = 0, as asserted. O 

Lemma B 2 . Let fA t B t c(2,y) be defined as in Eq. (26), 
and suppose there exists a nonzero element y* in GF( 8 ) 
such that 

fA l B ) c(x^y*) ~ 0, for all x £ GF( 8 ) (B- 2 ) 

Then, A = 0 and B = Cy\ Further, if Eq. (B-2) holds and 
if fA,B,c{Xiy) is not identically zero, then for any y ^ y* } 
the number of solutions x £ GF(8) to fA t B,c( x ,y) = 0 is 
exactly four. 

Proof: If Eq. (B- 2 ) holds, then by Eq. (26), one has 

Tr((v4y*)a; + (B + Cy*)x 6 ) = 0, for all x € GF( 8 ) 

(B-3) 

Then, since y* ^ 0, Lemma Bl implies that A = 0 and 
B+Cy* = 0, i.e., B = Cy* . This proves the first statement 
of the Lemma. To prove the second statement, assume 
that Eq. (B-2) holds and /a,b } c{ x > y) is not identically 
zero. Then, since it is already known that A — 0 and 
B = Cy* , it must be true that C ^ 0 , so that the equation 
/A,£,c(z>y) = 0 becomes 


Since C ^ 0 and y / y* , it follows that Cy * -f Cy ^ 
0, so that Eq. (B-4) has the form Tr (Dx 6 ) = 0, with 
D 7 ^ 0 . But since for z £ GF( 8 ), Tr (z) = 0 has ex- 
actly four solutions, viz., z = 0 , /?, /? 2 , /? 4 , it follows that 
Eq. (B-4) has exactly four solutions. □ 

Lemma B3. Let gD,E,F( x ,y) be defined as in Eq. (27), 
and suppose there exists a nonzero element y* in GF(8) 
such that 

9D t E t F{ x , y*) = 0, for all x £ GF(8) (B-5) 

Then, F = 0 and E = Dy * . Further, if Eq. (B-5) holds and 
if 9D,E t F( x > y) is not identically zero, then for any y / y* } 
the number of solutions x £ GF( 8 ) to gD,E t F( x ,y) — 0 is 
exactly four. 

Proof: The proof of Lemma B3 is similar to the proof 
of Lemma B2 and is omitted. O 

Lemma B4, Let A, J5, C, D, E> and F be elements 
of GF(8) such that 

fA,B,c{x,y) = gD,E,F(x,y) (B- 6 ) 

for all x G GF( 8 ), for two distinct values of y, say, y = y x 
and y = y 2 . Then, A = D, C = F, and B - E = 0. 

Proof: In view of the definitions in Eq. (26) and 
Eq. (27) of / and g, the given conditions are equivalent 
to 

Tr((Ay + Dy + E)x 

+ (B + Cy + Fy)x 6 ) = 0, for all x 6 GF( 8 ) 

(B-7) 

for y = yi,y2- Thus, according to Lemma Bl, Ay + Dy 
+ E = 0 and B -f Cy + Fy = 0 for y = y, , y 2 . The two 
equations Ayi+Dy { + E = 0 imply that A = D and E = 0, 
and the two equations Cy,- + Fy; +5 = 0 imply that C = F 
and 5 = 0 . □ 


Lemma B5. The code Aq defined in Eq. (28) is a 
TV ((Cy* + Cy)x 6 ) = 0 (B-4) [24, 12,8] code. 
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(B- 11) 


Proof: The mapping from 6-tuples [A, B } C } to, , ^ 2 ] 
to codewords in Aq is linear. The kernel of this mapping is 
the set of 6-tuples such that the corresponding codeword 
is 0, i.e. , 

fA,B,c(x,p r ) + e i =0, fori = 0,1,2 (B-8) 

for all x G GF{ 8). By substituting x = 0 into these equa- 
tions, one finds that c, = 0 for * = 0,1,2, so that in fact 

Ia,b,c{x,P 2 ) = 0, for i = 0, 1,2 (B-9) 

for all x € GF( 8). It then follows from Lemma B2 that 
A = D = C = 0. Thus, the kernel of the mapping contains 
only the 6-tuple [0, 0, 0, 0, 0, 0], and so the mapping is one- 
to-one. But since the set of 6-tuples is 12-dimensional, 
it follows that the code is also 12-dimensional. Thus, A 0 
is a (24,12) code. It remains to prove that its minimum 
distance is 8. 

To show that the minimum distance is 8, first consider 
the (24,9) code A' 0 defined by Eq. (28) with «o = «i = 
e 2 = 0. Each word in A' Q has three 8-bit segments, viz., 
the 8 bits corresponding to the function }a,b,c { x j J/) f° r 
y = P, P 2 , and p A . Since in each segment the bit corre- 
sponding to x = 0 is 0, each segment may in fact be viewed 
as a 7-bit codeword with components indexed by the con- 
secutive powers of a primitive root of GF{ 8). Thus, for the 
“ A,B,C ” codeword in A' 0 , the y-segment’s ith component 
is given by Ja,B,c{P' Aj)> where 

fA,B,c( x , y ) = Tr (Ayx + (B + Cy)x 6 ) 

= (Ay)x + (A 2 y 2 )x 2 + (jB 4 + CV)x 3 

+ {A A y A )x A + (B 2 + C 2 y 2 )x 5 + (B + Cy)x 6 

(B-10) 

It follows that each 7-bit segment is a codeword in the 
(7,6) binary cyclic code with generator polynomial g(x) = 
x — 1. In particular, each segment has even weight. The 
value of the weights modulo 4 can be computed by a 
theorem of McEliece-Solomon [9, Theorem 1] or [1, The- 
orem 16.33], which says that if an even-weight binary 
vector a = (ao, • ■ • , a »-i) * s described by its Mattson- 
Solomon (MS) polynomial (discrete Fourier transform) 
A(x) = A 0 + b v4 n _ix" -1 , i.e., if 


j=o 

where P is a primitive nth root of unity, and if ^(a) = 
E(=o 1)/2 AjA„-j, then w(a) = 2T 2 (a) (mod 4). The MS 
polynomials for the 7-bit segments are given by Eq. (B-10), 
and so the value of T 2 for the t/-segment is 

T a (y) = {Ay)(B + Cy) + ( A 2 y 2 )(B 2 + C 2 y 2 ) 

+ (B 4 + C 4 y A )(A 4 y*) 

= (AB + A*C*)y + ( A 2 B 2 + AC)y 2 
+ (A*B 4 + A 2 C 2 )y* 

= Tr((AB + A 4 C 4 )y) (B-12) 

If the three segments are combined into one 21-bit word, 
the overall weight is still even, and the overall weight mod 4 
is determined by the sum of the T 2 s, viz., 

r 2 = r 2 (/?) 4- r 2 (^ 2 ) + r 2 (/? 4 ) 

= Tr ((i4i^ -f A 4 c 4 )(l3 + $ 2 + p 4 )) 

= Tr(0) = 0 (B- 13) 

Thus, each 7-bit segment has weight 0, 2, 4, or 6, and the 
overall weight is divisible by four. Furthermore, if one of 
the segments has weight zero, then by Lemma B2 either 
the other two segments are both zero, or else the other two 
segments have weight 4. It follows that the weights in the 
(24,9) code are 0, 8, 12, and 16. Now the original code Aq 
is obtained from A f Q by complementing some or all of the 
segments, i.e., by replacing a segment of weight w with one 
of weight 8— w. Thus, in ^4o > segments have weight 0, 
2, 4, 6, or 8. But since 8 - w = w (mod 4), the weights 
in Ao must also be divisible by four, and so in Aq the only 
weights that can occur are 0, 4, 8, 12, 16, 20, and 24. The 
weight 4 can only occur as 0 + 0 + 4 and 0 + 2 + 2. Both 
of these cases can be eliminated by observing that since a 
zero- weight segment can only occur in an uncomplemented 
segment, and Lemma B2 says that if a codeword in Aq has 
a zero- weight segment, then either the other two segments 
both have weight 8, or both have weight 4. Weight 20 is 
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ruled out by observing that the complement of a word of 
weight 20 is a word of weight 4. d 

Lemma B5 says that the code Aq is a [24, 12,8] binary 
linear code. According to MacWilliams and Sloane [ 6 , Sec- 
tion 20.6], such a code must be equivalent to the Golay 
code. The next Lemma indicates that the code Bo defined 
in Eq. (29) is also equivalent to the Golay code. 

Lemma B 6 , The code B Q defined in Eq. (29) is a 
[24,12,8] code. 

Proof: The proof is virtually the same as the proof 
of Lemma B5. The key difference is that in place of 
Eq. (B-10), one has for the code Bq 

9 d,e,f(x, y ) = Tr ((Dy + E)x + ( Fy)x 6 ) 

= (Dy + E)x + (D 2 y 2 + E 2 )x 2 
+ (F 4 y 4 )x 3 + (D 4 y 4 + E 4 )x 4 
+ (F 2 y 2 )x 5 + ( Fy)x 6 (B-14) 

so that in place of Eq. (B-12) one has 

r 2 (y) = (Dy + E)(Fy) + (D 2 y 2 + E 2 )(F 2 y 2 ) 

+ (FV)(DV + E 4 ) 

= (EF + D 4 F 4 )y + ( E 2 F 2 + DF)y 2 

+ (E 4 F 4 + D 2 F 2 )y 4 

= Tr ((EF + D 4 F 4 )y) (B-15) 

and in place of Eq. (B-13) one has 

r 2 = W) + r 2 (p 2 ) + r 2 (/? 4 ) 

= Tr((EF + d 4 f 4 )(p + p 2 + 0 4 )) 

- Tr(0) = 0 (B-16) 

Further details are omitted here. O 


Lemma B7. The code A\ defined in Eq. (30) is a 
[24,9,8] code. 

Proof: Just as in the proof of Lemma B5, the mapping 
from 6 -tuples [ A , 0, C, eo, fi, 62 ] to codewords in A\ is a lin- 
ear, one-to-one mapping, which implies that A\ is a (24,9) 
code. Since A\ is a subcode of A 0 , its minimum distance 
must be > 8 . There are, however, many codewords in A\ 
of weight 8 , e.g., that obtained by taking A — C — 0 , 
eo = 1 , and = e 2 = 0 . □ 

Lemma B8. The code A 3 defined in Eq. (31) is a 
[24,5,12] code. 

Proof: Using the formula Eq. (26) for fA l B > c{x^y) } 
one finds that any codeword in A 3 can be represented as 
follows: 

[IY(C/?x 6 ) +f 0 |Tr(C7? 2 * 6 ) + e x \Tr(Cp 4 x 6 ) 

+ ( e 0 + t l)] x gG/’(g) (B-17) 

Two cases are considered: C — 0 and C ^ 0. If C = 0, 
then the codeword in Eq. (B-17) becomes [eo|ci|co 4* fi], 
which is either identically zero or has weight 16. If C ^ 0 , 
then since there are exactly 4 elements in GF( 8 ) with trace 
0, the codeword in Eq. (B-17) has weight 12. Thus, the 
only weights that occur in A 3 are 0, 12, and 16, and so A 3 
is a [24,5, 12] code, as asserted. O 

Lemma B9, The code A 4 defined in Eq. (32) is a 
[24,2,16] code. 

Proof: According to the definition in Eq. (32), each 
codeword in A\ has three 8 -bit segments. Either all three 
segments are identically zero, or else one segment is zero 
and the other two have weight 8 . Thus, in A\ the only 
weights that occur are 0 and 16, so that A 4 is a [24, 2, 16] 
code, as asserted. d 

Lemma BIO. Aq fl Bq = A\. 

Proof: Note that from Eq. (26) and Eq. (27), 

/a,o, c(s>y) = 9 A,o,c(x>y)- Thus, the code A 0 D B 0 con- 
tains the code Aj as defined in Eq. (30). To prove the 
opposite inclusion, note that by the definitions Eq. (28) 
and Eq. (29) of Ao and Bo, any word in the intersection 
will produce an equation of the form /a,s,c(^, P r ) + c, = 
gD,E,F(X)P 2 ) + for all x € GF( 8 ), for i = 0, 1,2. By 
substituting x = 0 on both sides of this equation, one gets 
€i = 6i, so that in fact, fA,B f c{x,P 2 ') = 9 D t E l F{^>P V ) for 
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all x 6 GF{ 8), for i = 0, 1,2. By Lemma B4, A = D, 
C = F, F = 0, and B = 0, so that a word in the inter- 
section must be of the form Eq. (30), i.e., it must lie in 


Lemma Bll. A\ H B\ = ^ 2 * 

Proof: Given the definitions in Eq. (30) and Eq. (33) 
of Ai and Bi, any word in the intersection j4i.fl Bi will 
produce an equation of the form fA,o,c( x t fl 2 ) + = 

9 o,e,f( x >P 2 ') + Si, for all x G GF(8) and t = 0,1,2. 
By substituting x = 0 on both sides of these equations, 
one gets e t * = 5,*, so that in fact one has fA,o,c( x ,P 2 ) - 
g 0 E F (x, /3 2 '). By Lemma B4, this implies A = E = 0 and 
C - F. Thus, the intersection j4 x fl B\ is exactly the same 
as A 2 , as defined in Eq. (34). O 

Lemma B12. A 2 0 B 2 = A 3 . 

Proof: Given Lemma B4 and the definitions Eq. (34) 
and Eq. (35) of A 2 and B 2 , this result is immediate. □ 


Lemma B13. A 2 fl B 2 — A 3. 

Proof: If a word in A 2 , as defined in Eq. (34), is the 
same as a word in B 2 , as defined in Eq. (39), then by 
setting x = 0, one finds that eo = <5o, ei = 5i» an d 6 2 = 

Thus, also /o f o,c(* > y) = 0o,E,a(*> v) for a11 * G Gir ( 8 ) 
and y = It then follows from Lemma B4 that 

C = E = 0, and so a word in the intersection ^2 H 
must be of the form described in Eq. (40). O 

Lemma B14. A 3 f) B 3 = A 4 . 

Proof: If a word in A 3 , as defined in Eq. (40), is the 
same as a word in B 3 , as defined in Eq. (41), then by 
setting x = 0, one finds that eo = Sq, c\ = $i } and c 2 — 
6 $ + Si. Thus, also /o, o t o(*,y) = yo.e.o^y) for all x G 
GF(8) and y = It then follows from Lemma B4 

that e = 0, and so a word in the intersection A 3 fl B 3 must 
be of the form described in Eq. (32). ^ 
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